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Abstract: 

In this paper, we study the detection of accelerator neutrinos produced at the China Spallation Neutron Source 
(CSNS). Using the code FLUKA, we have simulated the production of neutrinos in the proton beam on the tungsten 
target and obtained the yield efficiency, numerical flux, and average energy of different flavors of neutrinos. Fur¬ 
thermore, detections of these accelerator neutrinos are investigated in two reaction channels, the neutrino-electron 
reactions and neutrino-carbon reactions. The event numbers of different flavors of neutrinos have also been calculated. 
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1 Introduction 

In the past several decades, a number of spallation 
neutron sources have started to operate, such as the Los 
Alamos Meson Physics Facility (LAMPF), the Spalla¬ 
tion Neutron Source at Rutherford Appleton Laboratory 
(ISIS) [I], the Japan Accelerator Research Complex (J- 
PARC) [2|, the Spallation Neutron Source at Oak Ridge 
National Laboratory (SNS) etc. In recent years, 
new spallation neutron sources are under construction, 
such as the China Spallation Neutron Source (CSNS) 


0, the European Spallation Neutron Source (ESS) [§], 
etc. These spallation neutron sources are designed to 
provide wide multidisciplinary platforms for scientific re¬ 
search and industrial applications at national institutes, 
universities, and industrial laboratories @, [jj. The areas 
concerned include basic energy sciences, particle physics, 
and nuclear sciences. Table 1 shows the main technical 
parameters of several major spallation neutron sources 
at the GeV energy range. For J-PARC, the technical pa¬ 
rameters shown in the table refer to the Rapid Cycling 
Synchrotron (RCS). 


Table 1. Main technical parameters for several spallation neutron sources at the GeV energy range (ppp means 
protons per pulse). 



Extraction energy 
(GeV) 

Extraction power 
(MW) 

Repetition rate 
(Hz) 

Average beam Current 

H 

Intensity 

(10 13 ppp) 

Target 

LAMPF [8] 

0.8 

0.056 

120 

1 

2.3 

Various 

ISIS [9,101 

0.8 

0.16 

50 

0.2 

2.5 

Water cooled 







/Tantalum 

J-PARC 11.12] 

3.0 

1.0 

25 

0.333 

8.3 

Mercury 

SNS [13,14] 

1.0 

1.4 

60 

1.6 

16 

Mercury 

CSNS-I(II) [15,16] 

1.6 

0.1(0.5) 

25 

0.063(0.315) 

1.56(7.8) 

Tungsten 

ESS [17] 

2.0 

5.0 

14 

62.5 

110 

Tungsten 


The CSNS consists of an 80 MeV proton linac, a 1.6 
GeV RCS, a solid tungsten target station, and various in¬ 
struments for applications spallation neutrons fl8f The 
accelerator operates at 25 Hz repetition rate with an ini¬ 
tial design beam power of 100 kW and can be upgraded 
to 500 kW. As the exclusive spallation neutron source in 
developing countries, CSNS will be among the top four 
of such facilities in the world upon completion. Table 2 
shows the main design parameters of CSNS-I and CSNS- 



A large number of neutrinos can be produced at the 
beam stops of high intensity proton accelerators. They 
form the neutrino beams for basic scientific studies to 
better understand properties of neutrinos and better 
probes of the weak interaction force [Ii|. During the 
last several decades, many neutrino experiments have 
performed based on neutrino beams from spallation neu¬ 
tron sources j2o| and other proton accelerators. They in- 
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elude the Liquid Scintillator Neutrino Detector (LSND) 
fH at LAMPF, the Karlsruhe Rutherford Medium En¬ 
ergy Neutrino experiment (KARMEN) fill fill at ISIS, 
the Tokai-to-Kamioka experiment (T2K1 |23j at J-PARC 
and so on. For CSNS, similar to other accelerators, an 
intensive beam of accelerator neutrinos can be brought 
by the proton beam hitting on the tungsten target. By 
using a neutrino detector similar to MiniBooNE 24], dif¬ 
ferent flavors of neutrinos can be detected. Therefore, 
neutrino properties, such as neutrino mixing parameters 
and the mass hierarchy may be studied. 


Table 2. Main design parameters of CSNS. 


Parameter / unit 

CSNS-I 

CSNS-II 

Beam power on target/MW 

0.1 

0.5 

Linac energy/GeV 

0.08 

0.25 

Beam energy on target/GeV 

1.6 

1.6 

Average beam current///A 

62.5 

315 

Pulse repetition rate/Hz 

25 

25 

Ion type, source linac 

H~ 

H~ 

Protons per pulse/10 13 

1.56 

7.8 

Target material 

Tungsten 

Tungsten 

Target number 

i 

i 

Target size/mm 3 

50 x 150 x 400 

50 x 150 x 400 

Beam section/mm 2 

40 x 100 

40 x 100 


2 Accelerator Neutrino Beam 


The idea of using accelerator neutrino beam to study 
neutrinos was initiated independently by Schwartz and 
Pontecorvo, and the experiment was first carried out 
by Lederman, Schwartz, Steinberger and collaborators 
[19l I25| . Low energy neutrino beams can be produced 
by the decays of 7r and p at rest, which are generated 
in low energy proton beams hitting on targets. In a 
beam dump experiment, the target for the primary pro¬ 
ton beam, where the neutrino parent particles emerge, is 
the medium for absorbing or stopping the hadrons. Then 
no drift space is provided for hadrons to decay in. There¬ 
fore, high intensity and low energy proton accelerators 
(Pbeam ~ 1 GeV/ c) such as spallation neutron sources are 
commonly used [26l [271. The two recent experiments at 
a beam dump, LSND and KARMEN, have given contro¬ 
versial results on neutrino oscillations which have been 


resolved by the MiniBooNE experiment 28|, 

The dominant decay scheme that produces neutrinos 
from a stopped pion source is |30} 


• A* + + Vp. 


r™ = 26ns, 


( 1 ) 


The bulk of the 7r _ ’s generated are strong absorbed 
by the target before they are able to decay, and most of 
the fi ~’s produced from the n~ decay are captured from 
the atomic orbit, a process which does not give rise to 
v e . Therefore, the yield efficiency of n e is a factor of 10 -3 
to 10 -4 lower, and it will be neglected in the following 
discussions. 


Fragments 



Fig. 1. The production mechanism of accelerator 
neutrinos at CSNS. 


The CSNS beam stop will provide a copious flux of 
neutrinos, primarily from 7r + and /i + decays. Fig. 1 
shows the scheme of neutrino production in the tung¬ 
sten beam stop. By using the code FLUKA 31] and the 
main design parameters given in Table 2, the processes 
that the 1.6 GeV proton beam hitting the tungsten tar¬ 
get were simulated. The results show that, after the 
complete decays of 7r + and p + , three different species 
of neutrinos (i/ e , produced have the same yield 

efficiency which is about 0.17 per proton. In Table 3, 
the event numbers per year of accelerator neutrinos pro¬ 
duced at CSNS are given. We find that the three species 
of accelerator neutrinos have the same yield event num¬ 
ber per year which is 0.21xl0 22 for CSNS-I and 1.05xl0 22 
for CSNS-II. Since it will rise ultimately to 500 kW in 
the future with more neutrinos generated, the main pa¬ 
rameters for CSNS-II will be considered mainly in the 
following research. 


Table 3. Event numbers per year of accelerator 
neutrinos produced at CSNS. 



CSNS-I 

CSNS-II 

Proton number per year/10 22 

1.23 

6.15 

v & number per year/10 22 

0.21 

1.05 

v^ number per year/10 22 

0.21 

1.05 

v^i number per year/10 22 

0.21 

1.05 


The numerical flux of each particle (7r + , /x + , v e , v tl . 
v^) can be found in [321 ] 


followed by 

H + — >e + +V,j, + v e , T^ = 2.2ps, 
where r^(r M ) is the life time of 7r + (/i + ). 


( 2 ) 




ip(num/year) 
4:TrL 2 (cm 2 ) ’ 


(3) 


where L is the distance of the spallation target from the 
neutrino detector, and tp is the yield event number per 
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year. By using FLUKA, the yield efficiency of different 
particles can be obtained, then the numerical flux can be 
calculated using Eq. ©• Fig. 2 shows that: 

(i) Numerical fluxes of various particles involved de¬ 
crease with the distance L; 

(ii) u e and have the same numerical flux; 

(iii) The numerical flux of is larger than those of 
i/ e and Vp, 

(iv) The numerical flux of 7r + decreases very quickly 
with the distance L. 




L 


Fig. 2. The numerical flux of each particle (7r + , y + , 
v e , i/fj., as a function of the distance L. (a) n + 
and fi + -, (b) u e , v^, and 


From the simulation results of neutrino productions, 
the average energy of different species of neutrinos can 
be obtained. The energy spectra of the different speies 
of neutrinos are given in Fig. 3. It can be obtained that 
the average energy of v e is much smaller than that of 
and Dp. 



Fig. 3. Accelerator neutrino energy spectra at CSNS. 


Because of the flavor mixing there are oscillations 
among v e , and v r 33l- 351. In the three flavor mix¬ 
ing scheme, neglecting the matter effect, and using the 
fact Am 21 = 7.5 x 10 -5 eV 2 and L/E v sc 1, hence 


sin 2 (A , m| 1 -L/AE,,) sc 0, we can write the oscillation prob¬ 
abilities [36| 

P{y e — » Vy) ~sin 2 023 sin 2 20i 3 sin 2 (ATO 32 -L/4E„), 

P{y e — > v T ) cs cos 2 023 sin 2 20i 3 sin 2 (Am 32 -L/AE V ), (4) 
P(V M —>• v T ) ~ cos 4 0 13 sin 2 2023 sin 2 (A?n 32 • L/4E„), 


where | A to 32 |= 2.4x10 -3 eV 2 , sin 2 0 2 3 = 0.446, sin 2 0 13 = 
0.0237 [37|, E„ is the neutrino energy, and L is the dis¬ 
tance of the neutrino source from the detector. Given 
the distance L = 60 m the simulation results show that 
the average energy of i/ e , and are 33.0 MeV, 50.0 
MeV, 47.9 MeV respectively. The oscillation probabil¬ 
ities can be calculated as P(v e —> v^) sc 1.27 x 10 -6 , 
P{u e -»■ v T ) sc 1.57 x lO" 6 , and P(i/ M v T ) sc 1.26 x 10" 5 . 
Therefore, due to the very short baseline, the oscillations 
among ^ e , u fJ . and v r can be neglected and will not enter 
our discussions below. 

In the next section, processes of the accelerator neu¬ 
trino detection will be studied, and the corresponding 
neutrino event numbers observed through various reac¬ 
tion channels will be calculated. 


3 Detection of accelerator neutrinos 

The event numbers per year N t of accelerator neutri¬ 
nos observed through various reaction channels i can be 
calculated following jjjgj 

Ni = (j){v/year/cm 2 )-(Ji(cm 2 )-N T , (5) 

where </> is the neutrino numerical flux given in Eq. ©, 
(Ti is the cross section of the given neutrino reaction, and 
N t is the target number. 

For the accelerator neutrino detection at CSNS, a de¬ 
tector similar to MiniBooNE j|24| is adopted. It consists 
of a spherical 803 tons fiducial mass of mineral oil (CH 2 , 
density 0.845 g/cm 3 ) and has a fiducial radius of 6.1 m, 
occupying a volume of 950 m 3 . The total numbers of 
target protons, electrons, and 12 C are 

N^ = 6.90xl0 31 , N^ e) =2.76x 10 32 , 

A r y C) = 3.45 x 10 31 . 

According to the discussions given in the previous 
section, the yield efficiency of is a factor of 10 -3 to 
10~ 4 lower than that of the other neutrino species, and 
hence can be neglected. Then, the reaction channel of the 
inverse beta decay will be neglected. For the neutrino- 
proton elastic scattering, due to the Cerenkov energy 
threshold and quenching, only the high energy part of the 
neutrino spectra can be observed and it has always been 
taken to be too small in number to observe. In addition, 
due to the the proton structure, protons in the neutrino- 
proton elastic scattering are too difficult to identify and 
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have large systematic uncertainty (39|- Therefore, the 
neutrino-proton elastic scattering will be not considered 
in this paper. 

By using the said detector at CSNS, two reaction 
channels will be used to detect the different species of 
neutrinos (v e , i/^, and P M ): 


(1) Neutrino-electron reactions 

v e + e~ — >v e + e~ (CC and NC), 
v^-he —+ e (NC), 
v^-t-e (NC), 


where CC and NC stand, respectively, for the changed- 
current and neutral-current interactions, producing re¬ 
coil electrons with energy from zero up to the kinemat¬ 
ics maximum. The neutrino events observed through 
these reaction channels can be identified by the signal 
of the recoil electrons which are strong peaked along the 
neutrino direction [ 4 ( 3 , 41|. and this forward peaking is 
usually used for experiments to distinguish the electron 
elastic scattering from the neutrino reactions on nuclei. 


(2) Neutrino-carbon reactions 

For the neutrinos and 12 C system, there are one 
charged-current and three neutral-current reactions: 
Charged-current capture of v e : 

i/ e + 12 CV 2 JV + e-, E th = n.MMeV, 
12 N^ 12 C + e+ + v„. 


Neutral-current inelastic scattering of v e , 17 , and v 

v e + 12 C —> 12 C* + v e , E t f l = lh.llMeV, 

v^E 12 C —> 12 C*+ v' fl , E th = 15.11MeV, 

^+ 12 C^ 12 C-* + p;, E th = 15.11MeV, 
12 C *^12 C + 1 . 


The charged-current events have the delayed coincidence 
of a /3 decay following the interaction. The neutral- 
current events have a monoenergetic 7 ray at 15.11 MeV. 
Therefore, the charged-current and neutral-current reac¬ 
tions on carbon can be identified and observed by the 
neutrino detector (42|, 0 ] ■ 

The effective cross sections of the above two re¬ 
actions, the neutrino-electron reactions 0 , 0 and 
neutrino-carbon reactions 0 0, are given in Fig. 4. 

The neutrino event numbers per year can be calcu¬ 
lated using Eqs. m and ©. Fig. 5 shows the neutrino 
event numbers per year observed through two reaction 
channels, changing with the distance L. We found that: 

(i) The event numbers per year of different species 
of neutrinos all decrease with the distance L for both 
the neutrino-electron reactions and neutrino-carbon re¬ 
actions; 


(ii) The total event number of accelerator neutrinos 
observed through the neutrino-carbon channel is much 
larger than through the neutrino-electron reactions; 

(iii) For the neutrino-electron reactions, N Vt .(CC + 
NC)>N^NC)>N^(NCy, 

(iv) For the neutrino-carbon reactions, N v ( NC ) > 
N^{NC) > N„ e {NC)-, N^(NC) and N Ve (CC+NC) are 
both larger than N^^(NC). 

(a) Neutrino-electron reactions (b) Neutrino-carbon reactions 




Fig. 4. The effective cross sections as functions of 
the neutrino energy, (a) the neutrino-electron re¬ 
actions; (b) the neutrino-carbon reactions, a = 
/t,t, x = 



Fig. 5. The event numbers per year of accelera¬ 
tor neutrinos observed through various reaction 
channels as functions of the distance L. (a) 
the neutrino-electron reactions; (b) the neutrino- 
carbon reactions. 

Table 4. Neutrino event numbers per year ob¬ 
served through various reaction channels for L = 
60 m. The detector efficiency and beam-on effi¬ 
ciency are both taken to be 50%. 


Reaction 

N U( , (CC) + N„ e (NC) 

N^(NC) 

No u {NC) 

ue~ 

263 

123 

54 

u 12 C 

824 + 426 

2580 

1005 


In the future, if the distance L, the detector effi¬ 
ciency, and the beam-on efficiency are defined, the neu¬ 
trino event numbers per year observed through various 
reaction channels can be calculated accurately. For ex¬ 
ample, suppose the distance L = 60 m, the detector effi¬ 
ciency and beam-on efficiency both at 50%, the accurate 
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neutrino event numbers per year are given in Table 4. It 
is clear that there are a large number of accelerator neu¬ 
trinos which can be used for measuring neutrino cross 
sections. 

4 Summary and Discussion 

In this paper, the accelerator neutrino beam has been 
studied in detail. With the code FLUKA, processes of ac¬ 
celerator neutrino production at CSNS from the proton 
beam on the tungsten target have been investigated, and 
the yield efficiency, numerical flux, average energy of dif¬ 
ferent species of neutrinos have been obtained. We show 
that, after the complete decays of i r+ and p + , three kinds 
of accelerator neutrinos have the same yield efficiency 
which is about 0.17 per proton. Therefore, they have the 
same yield event number per year which is 0.21 xlO 22 for 
CSNS-I and 1.05 x 10 22 for CSNS-II. 

The detection of accelerator neutrinos through two 
reaction channels, the neutrino-electron reactions and 


neutrino-carbon reactions, has been studied, and the 
neutrino event numbers have been calculated. It is found 
that the total event number of accelerator neutrinos ob¬ 
served through the channel of the neutrino-carbon reac¬ 
tions is much larger than that of the neutrino-electron 
reactions. 

In our calculation, the detector efficiency and beam- 
on efficiency were not seriously considered. In the future, 
in the completion of the design of a neutrino detector at 
CSNS, the detector efficiency and beanr-on efficiency will 
be given. Then, the detector errors will need to be re¬ 
duced in the calculation of the neutrino event numbers. 
Furthermore, the statistical errors and systematic errors 
on neutrino fluxes and cross sections also need to be con¬ 
sidered for the design of neutrino experiment. 


The author would like to thank B.-L. Young, X.-H. 
Guo, S. Wang, and S.-J. Ding for helpful discussions 
and support. 
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